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Thermal losses from the hot end to the cold end of a Stirling cycle regenerator due to 
thermal dispersion through the regenerator matrix may significantly degrade the 
performance of the machine. Because of poor access to the void spaces within the porous 
medium, no direct measurements of thermal dispersion have been made and dispersion 
models have been derived indirectly. This is done by measuring the overall thermal 
performance of the regenerator and subtracting off the energy transfer caused by molecular 
conduction and advected enthalpy flows computed from volume-averaged fluid velocity and 
temperature. In the current program, a large-scale porous matrix consisting of stacked wire 
screens with a porosity of 90% is installed in a flow rig which is operated in a Reynolds 
number range that represents Stirling engine regenerator flow. Experiments are conducted 
to measure turbulent transport of momentum at the exit plane using hot-wires. The 
relationship of such turbulent transport terms to the thermal dispersion term in the 
volumetric-averaged energy equation for the regenerator matrix is developed and the 
measurements are used to determine cross-stream thermal dispersion. A dispersion model 
based upon the measurements is proposed and compared with models documented in the 
literature. 
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Nomenclature 

specific heat of fluid phase, J/kg K 
correlation coefficient 
characteristic length, m 
hydraulic diameter, m 
permeability of porous matrix, m 2 
thermal dispersion conductivity, W/mK 

thermal conductivity of fluid phase, W/mK 

the number of points selected to evaluate eddy viscosity 

Peclet number =ud /a f 

radial locations, m 
temperature, K 
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velocity, m/s 

thermal diffusivity of fluid phase, m 2 /s 
thermal dispersion eddy diffusivity, m 2 /s 
eddy diffusivity of momentum, , m 2 /s 
dispersion coefficient 
density of fluid phase, kg/m 3 
anisotropy factor 


Introduction 

P orous materials have been widely used for regeneration of engines, such as Stirling engines. Typically, the 
Stirling engine regenerators are made of wire screens or random fibers. The regenerator allows the ideal Stirling 
cycle machine to have a cycle efficiency equal to the Carnot cycle efficiency. Regenerator flow losses and heat 
transfer through the matrix contribute significantly to degrade the performance of the Stirling cycle, however. An 
important degradation mode is via the net energy flux through the regenerators directly to the coolers. 1 It was noted 
that a 1% decrease in energy storage in an engine regenerator matrix might necessitate a 4-5% increase in thermal 
energy input to a heater (i.e. a 4-5% decrease in efficiency), assuming the same indicated power. A major 
contribution to this energy transfer has become known as enhancement in axial heat conduction due to dispersion 
within the working gas (a component of which is due to eddy diffusion). It is a result of interaction of the flow with 
the regenerator matrix. 

In porous media, “dispersion” is described as “heat transfer due to hydrodynamic mixing of the interstitial fluid 
at the pore scale.” 2 It can be derived from the energy equation using volume averaging techniques and is defined as a 
correlation of the spatial fluctuations of temperature and velocity. For a flow with a high Peclet number, the thermal 
dispersion conductivity tensor, k d , can be modeled as 

k d = Xk f Pe (1) 

The dimensionless parameter, X , called the dispersion coefficient, depends on the structure of the porous medium, 
the flow structure and the flow direction. However, for a homogeneous isotropic porous matrix, the component of 
thermal dispersion due to eddy transport might be isotropic and a main contributor of the thermal dispersion in 
cross-stream direction. The eddy component of dispersion is the focus of the present paper. 

Extensive investigations have experimentally and theoretically determined the dispersion coefficient. A majority 
of the research has focused on a packed-particle bed since it is a simulation of many of porous media in nature. For 
this type of porous medium, Yagi and Kunni 3 determined the coefficient, X , equal to 0.09-0.1. Hsu and Cheng 4 
carried out experiments in cylindrical and annular tubes filled with glass spheres and derived the radial dispersion 
coefficient as a function of porosity. More recently, the radial dispersion coefficient was determined by Metzger et 
al. 5 to be about 0.03. For packed beds, the characteristic length, d , is the particle diameter. 

The porous structures of Stirling regenerators (consisting of wire screen or random fibers) differ from those for a 
packed bed. Moreover, the regenerators have a high porosity (80-95%) compared to packed beds (less than 40%). 
Therefore, none of the above models should be used directly in Stirling regenerators. Koch and Brady 6 analytically 
derived models for fibrous media. Different expressions were given to calculate the coefficients in different flow 
regimes and dor different directions. Nakayama and Kuwahara 7 numerically modeled convective heat transfer 
through a number of square rods as microscopic structures of a porous medium and gave correlations of dispersion 
in both longitudinal and transverse directions. Both models include the porosity as a parameter. Hunt and Tien 8 
measured heat transfer through a wall to porous media with porosities up to 97% and determined dipersion 

coefficients in that region to be about 0.025. They used the square-root of the permeability, <>[k , as the 
characteristic length. Few researchers have addressed modeling of dispersion in Stirling machines. Gedeon et al. 9 
extracted axial dispersion coefficient from experiments taken under oscillatory flow conditions. Results were 
computed from the the fluid peak value of bulk mean velocity and the fluid average temperature over a cycle and 
overall heat losses. They derived correlations of effective axial conductivity for various porous structures commonly 
used in Stirling regenerators 9 !0 . 

Kaviany 11 summarized different experimental methods to determine the thermal dispersion. The experimental 
methods frequently used (e.g. Ref. 5 and 8) require a solution of the macroscopic energy equation with the thermal 
dispersion coefficient included as a parameter. One can obtain the dispersion coefficient by matching the measured 
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temperatures or the measured heat transfer coefficients to those predicted by the solution. To the authors’ 
knowledge, no efforts have been taken to assess thermal dispersion by measuring flow velocity and temperature at 
the pore level, directly. There exist two challenges. One is that due to the complex geometry of typical porous 
media, no access is available for any currently existing probes. The other challenge is that even we do have access, 
the evaluation of spatial correlations of flow and temperature fields requires a number of probes within the porous 
medium so that the readings are taken simultaneously. It does not seem practical. The present paper is attempting to 
evaluate the eddy component of dispersion by direct measurements. The spatial variations of the velocity and 
temperature fields are described in terms of temporal variations with the aid of the following theoretical approach. 

Theoretical Approach 

A volume averaging technique has been commonly used in porous media studies. Any microscopic quantity, W, 
can be decomposed into the sum of the volumetric average (macroscopic) quantity <W> integrated over a 

representative elementary volume (R.E.V.), and a spatial variation, W , 

W=(w) + w ( 2 ) 

where W and W are functions of both time and space and <W> is a function of time, only. One can obtain one single 
volume-averaged energy equation by applying this technique to the energy equations for both solid and fluid phases 
while assuming thermal-equilibrium between the two phases, 4 ' 7 

P f c p { u) • V(r) = v • (k 0 v(T) - pjCpf (fa)) (3) 


The first term on the right hand side is the stagnant conduction by molecular diffusion. It depends on the fluid and 
solid phase thermophysical properties, the turtuosity of the structure and the porous medium geometry. The second 
term can be rewritten as 


- l>fC g ,(fi),k t -V(T ) 


(4) 


Here, the tensor, k d , is called thermal dispersion conductivity. Correspondingly, thermal dispersion diffusivity for 
axial transport is defined as 


£ (5, 

' PfCpf 

dx 

It can be seen that the thermal dispersion diffusivity, e d , is determined by the spatial average of the correlation 
between spatial variations of temperature and velocity. 

In order to obtain this correlation, one can choose a R.E.V., locate with a high spatial resolution a number of 
probes for temperature and velocity measurements within it, then take their readings simultaneously. Afterwards, 
one can calculate the spatial average quantities over the R.E.V. at any instant in time and the spatial variations at any 

instant in time and any particular location. The needed term <uT > is computed by spatially averaging uT . Note 
that it theoretically is a function of time since no time averaging has been done. This proposed experimental method 
cannot be realized. 

In contrast to the spatial average technique used in porous media is the evaluation of turbulence quantities taken 
in free and wall-bounded flows (but not porous media) in which a time-average is used to characterize the 
turbulence field. Any quantity can also be defined as the sum of the time-averaged value integrated over a time 
period sufficiently large relative to the characteristic time scale of the fluid flow situation, and a temporal fluctuation 
value. 
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w=w+w' 


(6) 


W' is a function of time and space and W is a function of space, only. 

These time-averaged quantities can be measured for an open flow by locating a probe at individual points in 
space and recording the readings continuously in time. This type of measurement is frequently taken. This raises the 
question, “Is it possible to use time-averaged values to describe the needed volumetric -averaged values?” Efforts 
have been taken to make a connection mathematically between volumetric-averaged quantities and time-averaged 
quantities. 12 Combining Eqs. (2) and (6) allows writing the volumetric average value, <W> , decomposed into the 
volumetric average of the time averages and temporal fluctuation terms as follows, 

(w) = (w + W' S j=(w S j+(W') (7) 

where (W^ may be a function of time but the temporal variation of (W') would be small if time and space scales 

for averaging are large enough relative to turbulent eddy passing time scales and turbulence length scales, 
respectively. The spatial variation of any quantity W can be written as, 

W =w-(w) = w+w'-{ffi + (w))~w-(w'j + w' (8) 

For the spatial variations at individual spatial points, it is assumed that time and space averaging is commutative. 
We note that (w) is an average of W over the representative elementary volume pore space. It could change with 

time but since it is a space -averaged value, it is expected to be essentially time-invariant. 

Applying Eq. (8) to velocity and temperature components of the fluid phase, then applying the volume-average 

technique to obtain the volumetric average of the product of u and T inside the porous medium, one can have 

(^uT^j = ((u - (w) -I - u )(T — (t} + T )j=<uT >-<u ><T > + <uT > + <uT > + <uT > (9) 

The first and second terms can be calculated from single-point experimental measurements since they do not require 

simultaneous measurements. In most applications, instead of <uT >, <uT > is of interest where the overbar 
represents an average over a period that is sufficiently long to capture all eddy scales of importance. If one can apply 

the time-average to Eq. (9), the third and forth terms go to zero. Then <uT > is reduced to the sum of three terms, 


(uT^=<uT >~<u ><T >+<uT > (10a) 

Similarly, for the cross-transport thermal energy transport and momentum transport: 


(vT^=<vT >-<v ><T >+<v T > (10b) 


(uv) — < uv > — < u >< v > +< u v > (10c) 

Since we are interested in the eddy transport component and believe that it is the main contributor to the cross- 
stream dispersion, we now focus on the cross-stream transport. As the first stage of our study, we now use the 
momentum transport results (Eq.(lOc)) to discuss thermal transport (Eq.(lOb)) by imposing Reynolds analogy, 
which states that momentum and thermal eddy diffusivities are nearly equal. 

Because the effect of eddy transport is office dominant in the cross-stream dispersion, the knowing the eddy 
transport term of the cross-stream dispersion is valuable. One can conclude that the first and second terms may 
cancel one another. The experimental support will be shown in the Appendix. The experiment is reduced to measure 
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< u v > within the porous matrix, or (u'vj , because the spatial and temporal averages are interchangeable. Then, 
Eq. (10c) allows saying, 



( 11 ) 


Recognizing the difficulty of measuring <wV> within the matrix, one might ask “Can the measurements be 
done within the turbulence structure of the flow just exiting the matrix instead of within the matrix?” Turbulent 
flows within porous media are different from free turbulent flows or turbulent flows far downstream of porous 
media. In the case of free turbulent flows, the eddies are formed as a break down of a shear layer. The eddies then 
strain one another according to the turbulence cascade theory. In porous media, the eddies are generated 
continuously due to separation off local solid surfaces. In some respects, this is similar to near-wall boundary layer 
turbulence where production of turbulence by inner-layer bursting continuously produces turbulent eddies. A probe 
located within the porous medium (or a boundary layer), not only sees the decaying eddies left over from upstream, 
but also the eddies produced locally that might have a large contribution to the eddy transport. Llittle research has 
been done on presenting a physical explanation of turbulent eddy transport in porous media. Masuoka and Takatsu 13 
provided the most detailed description of transport through porous media in terms of vortex transport. Vortices 
formed within porous media were categorized into two types: void vortices and pseudo vortices. They note that the 
pseudo vortices are the main contributors to thermal dispersion and play the role of long-distance (large-scale) 
momentum transport. We expect that they have a lifetime of at least several turnovers. 

Since the pseudo vortex has a large scale it might be measured immediately downstream of the porous medium 
instead of within the porous medium. This is similar to the assumption of advected “frozen turbulence” in which the 
turbulence field is unchanged as it exits the matrix and is advected downstream. One still can capture a major 
contribution of dispersion due to eddy transport. As the scale of the pseudo vortex might increase with the porosity, 
this statement becomes stronger as we approach porosity values seen in Stirling regenerators. 

Experimental Measurements 


Experimental Setup 

The porous medium in the present test is a simulation of a Stirling engine regenerator matrix operating under 
representative flow conditions. It consists of two hundred round layers of 6.3 mm x 6.3 mm square mesh (4 mesh 
per inch), stainless steel 304, welded screens. Each layer is a disk of 190 mm diameter. The wire diameter, d w of the 
screens is 0.81 mm. The screens are stacked in such a way that each layer is rotated 45° from the adjacent layer to 
simulate random wire arrangement of the real Stirling regenerator. The axial spacing of the screens is such that the 
regenerator has 90% porosity. The screen is punched out of welded screen sheets in such a way that the location of 
the center of the disc is random with respect to the center of the mesh in which it resides, adding to the randomness 
of the matrix. A more detailed description of the porous matrix is given in Ref. 14. A summary of important 
geometric parameters is given in Table 1. Under unidirectional, steady flow operation, the flow is introduced to the 
regenerator matrix through a duct of 2600 mm length and a diameter that matches the disk diameter of 190 mm. The 
flow in this duct is mixed so that the velocity distribution of the entry flow is flat. A triple-sensor, hot-film 
anemometer probe (TSI 1299BM-20) is located downstream of the regenerator matrix exit plane. The experimental 
setup is described in Fig. 1. Three velocity components, u , v and w , are measured, for various radial positions 
from -38 mm to 38mm (±1.5 inch) about 0.0, which is the radial center of the matrix. A sampling frequency of 

10kHz is used. By doing so, the turbulent shear stress, uv , can be evaluated. To create a radial velocity gradient, 
d u 

— - and a non-zero uv , a plate with a 12.5 mm (0.5 inch) diameter orifice in the middle is inserted into the 
dr 

regenerator at the 20 th layer of screen from the exit plane. As a result, local radial eddy diffusivity of momentum 
can be obtained as: 


“ ~ , du 


£ M ( r i) = ~UV / 


dr 


( 12 ) 
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where — ll V 


- du 

and vary with position, r. 

dr 



Triple-sensor probe 


Fig. 1 Experimental Setup for turbulent transport measurements 


Table 1 Geometry of the Porous Media 


Parameter 

Value 

Wire Diameter, mm 

0.8 

Mesh Size, PPI 

4x4 

Hydraulic Diameter, mm 

7.2 

Permeability, m 2 

1.86X1 0‘ 7 

Porosity 

90% 


Data Processing 

Figure 2 (a) shows a velocity profile at an axial location of 3.2 mm (0.125 inch) away from the exit plane of the 
regenerator. Due to the irregularity of the regenerator matrix, the probe might be located immediately behind a wire 
or immediately behind a pore. Therefore, the velocity distribution is highly variable in the radial direction, which 

makes it difficult to accurately evaluate the local, radial velocity gradients and the shear stress u v . To avoid this 
situation attributed to the close proximity to the mesh wires, the probe is moved further away from the exit plane. 
The velocity profiles from 3.2 mm (0.125 inch) to 50.8 mm (2.0 inch) are shown in Fig. 2. It can be seen that the 
radial distribution of velocity is becoming much smoother, ensuring a more accurate evaluation of the local radial 

velocity gradient and shear stress, u v . 

We might expect that eddy transport in the measurement zone may differ from that immediately behind the exit 
plane. With all the information at the various axial locations, we are tempted to assume that the turbulent structure 
does not change too much within a short distance downstream of the exit plane. To check the suitability of doing this 
we seek a relationship between the eddy diffusivity and the axial distance from the exit plane, as e M = e M ( jc ) . To 
eliminate the effect of randomness, only points with smooth gradients are chosen for evaluation of the averaged 
eddy diffusivity, (e M ) , at each axial location: 


( f A ir ‘ ),N 

Fj 


(13) 


The parameter N is the number of selected points for data processing, indicated in Fig. 2 as points with symbols 
in circles. As an example, the velocities at the axial location 50.8 mm (2 inches) away from the exit plane are plotted 
in Fig.2 i). In this case, twenty eight points are considered to be the points with smooth velocity gradients. Points 
with weak gradients were not processed. We find that several of the neighboring points to those eliminated (due to a 
weak gradient) show measured shear stress values that are of the wrong sign relative to the local mean shear, giving 
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a negative eddy diffusivity. We consider those points erroneous as well because, we suppose, we are not able to 
capture the relevant length scales of the Reynolds shear stress u'v ' . Points that are well inside of a continuous mean 
shear gradient region pass both of these tests and are processed. After these points have been eliminated, the number 
of selected points for evaluating the averaged values, N, is equal to 25. These are cycled points in Fig. 2 i). By 
similarly processing data at the other Fig. 2 planes, we can see that the number of points processed decreases 
significantly as the measurement plane becomes closer to the matrix exit plane. 

To derive a model, we use the Prandtl mixing length hypothesis to compute the eddy diffusivity as the product of 
a characteristic length (taken to be the hydraulic diameter) and a characteristic velocity (taken to be the in-matrix 
average velocity), 


£ M( r i)~^d h U (14a) 

Here, U values are the individual local velocity values at the selected points and d h is the hydraulic diameter of the 
regenerator matrix. An average diffusivity coefficient is defined as, 

(X)~YS e uWld h U.)IN (15) 

r t 

Figures 3 a) and b) indicate the values of < e M > and < A > , respectively, for various axial locations. 

The flow we have been measuring is that which is emerging from between the wires of the matrix. It has some 
characteristics of free jets. Measurements by others of eddy transport indicate that the eddy diffusivity remains 
constant with axial position for the near-field of axisymmetric free jets. If the present case could be modeled as an 
array of free jets, we may be inclined to extract from Fig. 3 a) a constant eddy diffusivity of around 1.7xl0' 4 m 2 /s 2 
and from Fig. 3 b) a corresponding coefficient of 0.02. We put low emphasis on the first three date points because so 
few of the measurements pass the two acceptance criteria described above. 
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Fig. 2 Velocity Distribution at various axial locations 
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Fie. 2 Velocity Distribution at various axial locations 
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Fig. 3 Values Calculated by Eqs. (13), (15), (16) and (17) 

We next test the hypothesis that, although the turbulence levels may change from within the matrix to the exit 
plane and beyond, the structure does not. To test this assumption with our measured data, we define a correlation 
coefficient and an anisotropy factor as follows 


< mV >= C t 



(16) 


a = 


/ v 2 /Vm 2 


(17) 


If the flow structure does not change from within the matrix to the outside of the matrix, and in the immediately 
vicinity beyond, we might expect constantjvalues for these two. If so, the experiment can be reduced to measuring 

the correlation coefficient, C t , and the anisotropy factor, a , in the exit stream, and only the intensity value 
within the matrix or at the exit plane of the matrix. 


< u'V >=C t a < u 2 > (18) 

The correlation coefficient, C h and a, are computed for each x-plane. Radius-average values at each axial 
position using the selected points are substituted for volume averages in doing this evaluation (See Fig 3 c) and d)). 
This C t (x) is extrapolated to the x=0 position to give about 0.2, again with reduced emphasis on the first three data 
points. It is observed that a value of 0.9 for a is maintained along the streamwise direction. This supports our 
previous assumption that the turbulence structure does not change in this 50 mm length (again, advected “frozen 
turbulence”). 
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Discuss and Comparison 

From the current experimental measurements, a model is suggested for evaluating the eddy transport component 
of the cross-stream dispersion in stacked wire screens or similar porous structures. 

=0.02 d h U (14b) 

A comparison of the models by different researchers is given in Table 1. Since these models employ different forms 
and different length scales, for the purpose of comparison, all dispersion coefficients are computed based upon the 
hydraulic diameter of the current porous medium using the parameters listed in Table 1 for conversion and local in- 
pore velocity. 

It can be seen that the thermal dispersion coefficients scatter widely, even for the similar porous materials. We 
note, however, that the modeled conditions differ from case to case, as will now be discussed. 

The value from the direct measurements (Eq. (14b)), is computed from the last term of Eq. (10c) and it is 
assumed that the first two terms cancel. This assumption is supported by our measurements. If this assumption is 
not so, Eq. (14b) would be in error. Also, if the probe were too small to resolve the important scales of the turbulent 

flow, then mV and A would be low. Studies of the scales the probe might resolve indicate that scales on the pore 
size should be resolved by our probe but turbulence scales on the wire size might not be resolved. These scale 
resolution studies were conducted in a fully-developed turbulent pipe flow in which the turbulence structure is well 
characterized. It appears that the larger eddies responsible for dispersion are being resolved so that filtering due to 
probe size is not expected to be a major problem. 

Theoretical model given by Koch and Brady 6 is derived from an isotropic packed fibrous media which is the best 
presentation of the porous matrix in the present study. The model by Nakayam and Kumahara 7 is developed from 
numerical simulations using a microscopic structure of square rods. All the rods in their structure are in-line and 
perpendicular to the flow direction, which is similar to the aligned fibers studied by Koch and Brady 6 . Such a 
regular pattern might result in a reduced dispersion. 

Although the model given in Sage Manual 10 for predicting dispersion in Stirling regenerator is for a similar 
structure of woven screens, it is derived experimentally from an overall measurement of the axial heat conduction 
losses 9 . It should be used for the evaluation of streamwise dispersion. The values from the study by Hunt and Tien 
relate to eddy dispersion near a heated wall. We know that a wall tends to suppress turbulent mixing in boundary 
layers so we might expect that dispersion in a near-wall region may be reduced from dispersion values which apply 
out in the matrix far from the wall. This may explain the smaller value than given by . 


Table 2 A comparison of the thermal dispersion coefficients 


Researchers 

Models 

Porous Media 

Current researchers 

Sfof = 0.02 dfoU 

Welded Screen 

Koch and Brady 6 

= 0.02 d h U 

Isotropic Fibers 

Koch and Brady 6 

e M = 0.0004 d h U 

Aligned Fibers 

Nakayama and Kuwahara 7 

e M = 0.0018 

Square Rods 

Sage Manual 10 

e H =0.061 d h U 

Woven Screen 

Hunt and Tien 8 

e H = 0.0015 d h U 

Fibrous Media 


Conclusions 

An effort has been made to find a measurement technique for directly evaluating the eddy component of thermal 
dispersion. A relationship between thermal dispersion within a porous medium and turbulent transport in the exit 
flow has been discussed. The former cannot be measured with the current techniques while the latter can be 
measured in the laboratory. As a first step, eddy transport downstream of the regenerator matrix is measured at 
various axial locations. A model for the eddy component of cross-stream dispersion is derived from the experimental 
results. This model is compared to some other models for evaluating thermal dispersion or turbulent transport in 
porous structures. The comparison shows that the correlations scatter widely. However, we presented reasons why 
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some might not agree with others. More work is needed to directly evaluate the streamwise dispersion and eddy 
transport of thermal energy. 
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Appendix 

The last term of Eq. (10c) is measured directly. It is assumed that the first and second terms might cancel one 
another or that the sum of the two is relatively small compared to the last term in Eq. (10c). A comparison of the 
terms computed from the measured results at various axial locations is made in this appendix. 

In order to show the relative importance between the sum of the first two terms <uv > - <u ><v > and 

< u v > , the ratio of the two is plotted in Fig. A-l. As a reminder, N data points passed the two criteria discussed 
above and were thus processed to evaluate the eddy diffusivity. For each of these N points, one neighbor on each 
side plus itself are taken as a representative pore size to compute the local spatial averaged values 

<u > ,<v > ,<Iiv > and <u v > . It is noted that the points that are chosen as neighboring points might fall within 
the selected N points, or might not meet the selection criteria. In a case where all three points of the set meet the 
criteria, we declare it to be a good data point. As an example, for the plane x=50.8 mm, we have 15 good data points 

out of the 25 points. Plotted is the x-plane-averaged value of (<uv>-<u><v>)/<u'v'> computed from the 
data and averaged over all the selected N points. 
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The ratio seems to be nearly zero for the x-stations where there is a sufficiently large number of good data points 
(further from the exit plane of the matrix). For the first three axial locations we have few good data points and the 
values scatter widel. For the next three, the scatter is also large and the values are much nearer zero. The last three 
have many more good data points than those of smaller x-value, scatter very little and are nearly zero in value. We 
feel that they are the most reliable data points for this evaluation. Therefore, based on these current experimental 

data, we can argue that the sum of the first two terms in Eq. (10c) is small relative to the term(wv) and that the 
measured term < u v > is approximately equal to the term (uv ) . 
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Figure A-l. (<uv>-<uxv>)/<u'v> at various axial locations 
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